%% Comparison of Errors: 
JLF = LJMM_F;
JLS = LJMM_SC;

err = reshape(U_vec./u_vec - 1,[N,T]);errT = err(:,T);errT(isnan(errT)) =[];
err_FP = reshape(U_vec_FP./u_vec - 1,[N,T]);err_FPT = err_FP(:,T);err_FPT(isnan(err_FPT)) = [];
errJLF = reshape(JLF./u_vec - 1,[N,T]);
errJLS = reshape(JLS./u_vec - 1,[N,T]);


BBK = max(abs(( errT )));
BBK_FP = max(abs((err_FPT)));
LJF = max(abs((errJLF(:,T))));
LJS = max(abs((errJLS(:,T))));

BBK2 = norm((errT))/sqrt(size(errT,1));
BBK_FP2 = norm((err_FPT))/sqrt(size(err_FPT,1));
LJF2 = norm((errJLF(:,T)),2)/sqrt(N);
LJS2 = norm((errJLS(:,T)),2)/sqrt(N);

ErrorJL = [LJF, LJS,LJF2, LJS2];
ErrorBBK  = [ BBK, BBK_FP, BBK2, BBK_FP2];